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The Hall viscosity describes a non-dissipative response to strain in systems with broken time- 
reversal symmetry. We develop a new method for computing the Hall viscosity of lattice systems 
in strong magnetic fields based on momentum transport, which we compare to the method of 
momentum polarization used by Tu et al. [Phys. Rev. B 88 195412 (2013)] and Zaletel et al. [Phys. 
Rev. Lett. 110 236801 (2013)] for non-interacting systems. We compare the Hall viscosity of square- 
lattice tight-binding models in magnetic field to the continuum integer quantum Hall effect (IQHE) 
showing agreement when the magnetic length is much larger than the lattice constant, but deviation 
as the magnetic field strength increases. We also relate the Hall viscosity of relativistic electrons in 
magnetic field (the Dirac IQHE) to the conventional IQHE. The Hall viscosity of the lattice Dirac 
model in magnetic field agrees with the continuum Dirac Hall viscosity when the magnetic length 
is much larger than the lattice constant. We also show that the Hall viscosity of the lattice model 
deviates further from the continuum model if the C 4 symmetry of the square lattice is broken to 
C 2 , but the deviation is again minimized as the magnetic length increases. 


I. INTRODUCTION 

The topological response properties of the quantum 
Hall effect have been intensely studied for more than 
three decades, and begun with the understanding that 
the quantized integer/fractional Hall conductance itself is 
a topological phenomenon.The field has since under¬ 
stood that quantum Hall systems (with and without mag¬ 
netic helds) also exhibit remarkable responses to changes 
in geometry. 

One interesting piece of the geometric response is the 
non-dissipative Hall viscosity.^ The Hall viscosity r]H 
is an off-diagonal response coefficient that is only non¬ 
vanishing when time-reversal symmetry is broken. Un¬ 
der time-dependent shear strain, the viscosity tensor ry 
relates the stress tensor T to the strain rate it: 


formed via the adiabatic transport of the Hall fluid un¬ 
der shear strain on a torus.For Schrodinger elec¬ 
trons at integer filling factors, this type of calculation 
yields tjh = hvpjA where p is the electron number den¬ 
sity, and V is the integer filling fraction.More re¬ 
cently, Ref. 17 developed Kubo formulas for the Hall 
viscosity which obtain the same result. Also, a new pos¬ 
sibility for calculating the Hall viscosity was proposed 
via the so-called momentum polarization entanglement 
technique,though there is very little explicit discus¬ 
sion of the results of this method in the literature (see 
Ref. 39 for a very recent article). Remarkably, from the 
adiabatic transport calculations it has been shown that 
for rotationally-invariant integer and fractional quantum 
Hall systems in large magnetic fields, the viscosity is 
quantized in units of the density^*^ and takes the form 
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where the strain tensor is constructed from a sym¬ 
metrized gradient of the local displacement If only 
dissipative viscosity coefficients are present, e.g. the bulk 
and shear viscosities, then and is thus 

symmetric under exchange of {pv) with {ajf). However, 
being non-dissipative, the Hall viscosity generates an an¬ 
tisymmetric piece satisfying 


tiiyoiB ctBni/ 

Oh = -Oh ■ 


When the system is 2D and isotropic, the antisymmet¬ 
ric part of the viscosity tensor is determined by a single 
parameter pn that gives, in an isotropic 2D orthonormal 
frame:^’® 


„1112 _ „1222 
Oh — Oh 


Vh 


Oh^^ = 0 . 


The Hall viscosity can be calculated using a variety 
of different methods. The first calculations were per¬ 


OH = '^hp (1) 

where k is a universal number characterizing the partic¬ 
ular integer/fractional quantum Hall phase, and p is the 
uniform electron number density. Generically, the Hall 
viscosity has units of [^] for some length scale £, but it 
need not always retain such a clear quantization in terms 
of the particle density. 

The goal of this article is two-fold: (i) we introduce 
a new method for the calculation of the Hall viscos¬ 
ity using momentum transport, and compare with the 
extraction of the Hall viscosity via the momentum po¬ 
larization entanglement method, and (ii) we study the 
properties of the viscosity in two different lattice realiza¬ 
tions of the Landau-level integer quantum Hall problem 
(square-lattice Hofstadter, lattice-Dirac model in a mag¬ 
netic field), and illustrate the competition between con¬ 
tributions of the viscosity from the lattice-length scale 
and the magnetic-length scale. Our article is organized 
as follows. In Section H, we present the momentum trans¬ 
port method used here to compute the Hall viscosity. We 
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also review the momentum polarization method, which 
allows computation of the Hall viscosity from the entan¬ 
glement spectrum. In Section III, we describe the appli¬ 
cation of both methods to the continuum Landau level 
problem, and compare the results to previous calculations 
of the Hall viscosity. In Section IV, we present numerical 
calculations of the Hall viscosity for the Landau levels 
of a tight-binding model (the Hofstadter model) and dis¬ 
cuss the results. In Section V we calculate the Hall vis¬ 
cosity of the Landau levels of the continuum Dirac equa¬ 
tion, with and without a mass term, building on previous 
work by Kimura^®. Finally, in Section VI we use a lat¬ 
tice analog of the continuum Dirac system in a magnetic 
field, and study the Hall viscosity for comparison with 
the continuum results. 


II. METHODS 

We will consider two independent methods for calcu¬ 
lating the Hall viscosity in our example systems. The 
first method considers the transverse flow of momentum 
when a cylinder is strained with an area-preserving de¬ 
formation. The second method uses the entanglement 
spectrum to calculate the phase acquired by the many- 
body wavefunction when half of a cylinder is sheared; 
from this phase, one can extract the central charge^® 
and Hall viscosity.^^ Let us introduce and review both 
of these methods. 

a. Momentum Transport For the first method we 
study the off-diagonal components of the stress tensor, 
which represent the momentum flux. We will write the 
geometric deformation in terms of the strain tensor. To 
be explicit, if Ua is the displacement vector, then, to low¬ 
est order, the strain is^° 

'Uq ,/5 — T • 

In terms of the geometry, if ds and ds' are the original and 
deformed length elements, respectively, then to lowest 
order in the deformations, 

ds'^ = ds^ + 2uapdx°‘dx^. 

Due to the structure of the Hall viscosity terms in the 
viscosity tensor, shear strain causes momentum transport 
in the direction of momentum, but pressure/stretching 
causes momentum transport orthogonal to the direction 
of momentum. We find the latter is more easily studied in 
lattice systems when using a cylinder geometry, therefore, 
we consider deformed metrics of the form 

ds^ = ^dx"^ + a^dy^ ( 2 ) 

where a can vary. This deformation is area-preserving 
(shear), so we need not isolate our momentum transport 
results from effects induced purely by changes to the den¬ 
sity. 


To calculate the momentum transport under this de¬ 
formation let us consider a cylinder which is periodic in 
the j/-direction with a circumference Ly. As a varies, the 
strain rate is 

• ^72 *72 

u = - 7zdx -\-aady . 

a'^ 

The corresponding stress tensor components arising from 
the Hall viscosity contributions are 


T^y = Ty^ = -2r]H-- 
a 

Note that when the system is anisotropic (a ^ 1) we can 
have 7 ^ 0; unfortunately that term will not appear 
in this component of the stress tensor, so it cannot be 
extracted; but hence, it also will not affect our calculation 
of the other viscosity coefficients. 

As the metric is deformed we want to study the amount 
of momentum transported from the left-half of the cylin¬ 
der to the right-half. Consider cutting the cylinder at 
X = Xcut; if Tr is the projection operator onto the right 
side of the cut, then the total momentum on the right of 
the cut is 


{Py^R) 




dy Hj^ 


where Uy is the momentum density. Typically, we will 
choose cCcut = 0 with the cylinder placed symmetrically 
around this point. The stress tensor gives the momentum 
flux across the cut, i.e. 


^ {PyPR) = - ^ dyTy-= 2LyyHaa. 

From this equation we can immediately read-off the im¬ 
portant result: 


Vh 


1 d 

Ly da^ 


{PyPR) ■ 


( 3 ) 


We will use this relationship between rjH and the a- 
dependence of the half-cylinder momentum to calculate 
the viscosity. 

b. Momentum Polarization The second method we 
discuss uses the entanglement spectrum^^d2 determine 
the momentum polarization.^^ The momentum polariza¬ 
tion was initially proposed to calculate the topological 
spin and central charge of the conformal field theory at 
the edge of a topological phase. For a system in a cylin¬ 
drical geometry, these data are extracted from the expec¬ 
tation value of the operator Ty , which globally translates 
the left half of the cylinder in the periodic direction. The 
expectation value can be computed using the reduced 
density matrix,®® 

X={G\T^\G)=irR{pLT^) (4) 

where [G) is the ground state. Ref. 38 shows that A can 
be easily calculated for free-fermion systems using the 
entanglement spectrum. 
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To see that the topological spin and central charge can 
be extracted from this expectation value, consider that, 
in the long-wavelength limit, the reduced density matrix 
of a cylinder cut in half can be written in terms of the 
Hamiltonians Hu and of the respective conformal 
edge theories of the left and right edges of the left half¬ 
cylinder only:^® 

PL = PLl ® PLr = 


The relevant half-cylinder translation operator is 


Ty = exp 


^ {Pi + Pr) Ay 


cylinder, but we will find it more useful to compute this 
formula in terms of the projections onto the right half, 
Cn,ky = 1 — • Using these identities, it is convenient 

to rewrite (5) as 

^= n ^ [(1 ++ (1 - ( 2 ^'=-" - 1 )] • 

n,fcy 

( 7 ) 

In a remarkable extension of this work. Ref. 22 shows 
that for quantum Hall states one can extract the Hall 
viscosity from the imaginary part of the “non-universal” 
coefficient a. Explicitly they find 


where Ay is the distance translated (which we take to be 
a multiple of the lattice constant for lattice systems), and 
Pi and Pr are the generators of translations (momentum 
operators) of the left and right edge theories on the half¬ 
cylinder, respectively.^®’^^ Since the left-most edge is far 
from the right half, /3/ —> oo, and only the ground state of 
the left edge contributes. The ground state expectation 
value of Pi is /i—c/24 where h is the topological spin and c 
is the chiral central charge.Therefore, the contribution 
of the left edge is^® 


/ 

27ri ^ 

\ 

27rf . c \ 

1 PLl exp 

-j—PiAy 

1 = exp 



On the other hand, Pr takes a finite value because the 
right edge is entangled with the right half-cylinder. In 
general, the right edge gives a non-universal contribu¬ 
tion®® 


ULr 


(^PLr exp 


27ri 

. T'y 


PrAy 


exp [—Lya]. 


From this we see that one can extract the central charge 
and topological spin. 

For free fermions, A is easily calculated in terms of the 
entanglement spectrum for a cylinder by the formula®® 


^ = n b -f (1 - tanh 

n,ky 


( 5 ) 


where ]/[„ is a product over the bands and y-momenta, 
and ^ky,n is the entanglement eigenvalue of the state in 
band n with momentum hky. The entanglement eigen¬ 
values can be expressed in terms of the eigenvalues of the 
free-electron, equal-time correlation function, 

1 - 

a„,n=l0g (6) 

'^ky ,n 

where are the eigenvalues of 

where ky are the momenta in the periodic direction, i,j 
run-over the lattice sites on the left half of the cylinder, 
and a, b run-over all of the onsite degrees of freedom. 
Note that this projects states onto the left half of the 


A = exp 





( 8 ) 


where additional non-universal terms that scale differ¬ 
ently with Ly have been dropped. In their work they 
consider a full twist such that Ay = Ly, but the result 
carries over for smaller Ay as well. Thus, the viscosity 
and central charge can be extracted from a fit of LyArgA; 
the former from the quadratic coefficient, the latter from 
the constant coefficient. 

We can understand how the momentum polarization 
phase encodes the viscosity by considering the action of 
the shear strain generators on the ground state. Here, 
we will show that the Hall viscosity can be extracted by 
comparing the momentum polarization calculated with 
a real-space cut to the phase taken with an orbital cut 
following Ref. 34. We note that Ref. 34 identified two 
distinct contributions to the Hall viscosity, and the con¬ 
tribution which interests us here is due to changing the 
shape of the Landau orbitals under shear strain. The sec¬ 
ond contribution, the guiding center Hall viscosity, comes 
from the electron correlations and is absent in the inte¬ 
ger quantum Hall models we study here. We will review 
how the momentum polarization phase calculated with a 
real-space cut encodes both Hall viscosity contributions. 
Although we consider only the integer effect, the guid¬ 
ing center Hall viscosity also has a super-extensive term 
due to the non-zero net momentum in each half of the 
system.®"* We will show that this background can be sub¬ 
tracted by calculating the momentum polarization phase 
with an orbital cut and comparing the two results. 

For most of the remainder of this section we closely 
follow Ref. 34. First, let us decompose our physical 
coordinate R into a guiding center coordinate r and an 
orbital coordinate r: 


R = r -I- r. 


There is a metric associated with the physical coor¬ 
dinate R, as well as metrics and associated with 
each coordinate r and r, respectively. Let the operators 
generate shear strain (area-preserving deformations) 
associated with the metric likewise, let be the 
shear strain generators associated with 5 ^^. These gen- 
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erators obey commutation relations^"^ 


Now, using the strain field in Eq. (10), we find 


[A^^ A“^] 
V'', A“^ 
A“^ 


= ^ (e^“A"^ + ^ v) 

= 0 . 


( 9 ) 


The strain generator in the physical coordinate is 


Ares = ('I'l U{a) |4') 

= (4'| u{a) |4') (Tl u{a) |4'), 

where Ares is the momentum polarization phase A in 
Eq. (7) computed with the real-space entanglement spec¬ 
trum. Now, the expectation value of u is the momentum 
polarization phase computed with the orbital entangle¬ 
ment spectrum.while the expectation value of u is 


so that the unitary operator implementing strain on 
quantum states is^’^’^^ 


U(a) = exp 




where a^i, is a symmetric matrix parametrizing the 
strain. Because the strain generators on each coordinate 
commute, we can also write this as the product of strain 
transformations on each coordinate: 


U(a) = u{a)u{a) 


u{a) = exp 
u{a) = exp 


i J fi^Rap^(R)A‘ 

i f ci^Rap^(R)A‘ 


To first order in the variation in the metric under 
strain is^^ 


('I'l u{a) |4>) = (4>| exp 
= (4>| exp 


= (4>| exp 


i J (i^RQ;^^(R)A'' 

z J d^Ra^^i-R)X^ 
i J d^R6{x)AyX^ 


IVI/) 
Ivf-) 


|vl/) 


= exp 


i / d2Rj(a:)Ay(4-|A“^|4') 


= exp[iLyAj/(^|A"n^)] 


where we have kept terms only to first order in Ay. Using 
the strain generator commutation relations in Eq. (9), we 
substitute 

(4-1 A^^ 14-) = i (4-1 [A^^, A^^] 14-) = = -^yH 

to find 


'^G'/xi^(R-) — (R)Q:/3!/(R) -f // O t/. 

In our particular case, where we shear half the cylinder, 
this gives 

y)= 0 ^ 


The momentum polarization expectation value A (c.f. 
Eq. (7)) is just the ground state expectation value {U{a)) 
under this strain field. 

Before we proceed to compute the required expectation 
values and find the momentum polarization phase, let us 
see how the Hall viscosity enters the calculation. We can 
represent the viscosity tensor in terms of the adiabatic 
curvature of the ground state under shear strain:^’^"^ 


(R) = 27 Im 


(ivI'(Q;) 

dcXfjii/ (H') 


(74-(a) 
do.oij3 (I^) 


= -ih (4-1 [A^'^, A“^] 14-). 


where Ivk)®)) = U{a) Ivk). Because the strain generator A 
is the sum of orbital and guiding center strain generators, 
we conclude that the viscosity also has contributions due 
to each strain generator, which we separately denote 


^ M -^ a / 3 ( R ) 

~Mi.a/3(R) 


—ih (4/1 
—ih (4-1 


A“^] 14-) 

A“^l 14-). 


(d/j u{a) 14/) = exp 


--LyAyyn 


Returning to our expression for the momentum polar¬ 
ization phase, we have 


Ares = Aoes exp 


--LyAyyn 


where Aqes = (4/| m(q;) |4/) is the momentum polariza¬ 
tion phase computed with the orbital entanglement spec¬ 
trum. Hence, we can determine that an alternate form 
of the (orbital contribution to the) Hall viscosity is given 

by 


Vh = - 



Ares 

Aoes 


( 11 ) 


for systems in uniform magnetic fields. 

Now that we have introduced the two separate meth¬ 
ods for calculating the viscosity we will apply them to two 
different continuum systems, and their matching lattice 
regularized models. 


III. CONTINUUM LANDAU LEVELS 

Let us begin with the conventional Landau level prob¬ 
lem of 2D electrons in a uniform magnetic field, and 
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FIG. 1. The Hall viscosity (tih) of the specified integer quan¬ 
tum Hall states calculated by the momentum transport (up¬ 
per panel) and momentum polarization (lower panel) meth¬ 
ods. The calculation converges when Lx > 2 ^/^£b, i.e. when 
each half of the cylinder is wider than a single wavefunction. 
The Hall viscosity is given in units of fipo where po — II2 tt£% 
is the electron density of the lowest Landau level. 


consider the possibility of geometric deformations sim¬ 
ilar to Ref. 7. The Hamiltonian for electrons in a 
background electromagnetic field subject to the metric 
of Equation (2) is 

On a cylinder which is periodic in y (with circumference 
Ly), and with a uniform magnetic field normal to the 
cylinder (in the Landau gauge, = 0 and Ay = Bx), 

As is conventional, we define the lowering operator 


, : ' 
\/2heB 


aPx - {Py + eBx) 


a 


( 12 ) 


and its adjoint, ab It is easy to verify that their commu¬ 
tator is 


so that these are the usual ladder operators of quantum 
harmonic oscillator. The Hamiltonian is 


H = huj 


a^a - 


where co = eB/m is the cyclotron frequency. 

The lowest Landau level wavefunction satisfies 


a(o) 


= 0 


a<l>k,c 

where we are using py = hk. The raising operator 
generates the higher Landau levels, 

1 




(at)' 


e- 


The general formula for the wavefunctions for the n-th 
Landau level is 


exp 


4 "a( 3 ^. y) = 


ikv- 


(2^n\a£BLyy/T^)' 


Hr, 


kil 


ais 


with k = 2Trn/Ly for n G Z 


(13) 

and where the magnetic 
length Hn is the n-th Hermite polynomial. 

When a = 1, i.e. in the absence of any metric de¬ 
formation, the wavefunctions assume their well-known 
isotropic form. 

Let us now present the calculations for the Hall viscos¬ 
ity using the two methods we presented in the previous 
section. To calculate the Hall viscosity by the momentum 
transport method at a filling n, we need only compute the 
derivative (with respect to a^, c.f. Eq. (3)) of 


v-l K 


{Py'PR) = Y. E 


(") 


n^Ok^-K 


poo pL 

where = dx 

Jo Jo 


dy 


(a:, y) 


(14) 


(15) 


We note two things: (i) is just the probability of 
finding a particle on the right {x > 0) half of the cylin¬ 
der, given that the particle is in the state ^[,4) and (ii) 
these quantities match the correlation-function eigenval¬ 
ues Ck,n if one calculates the entanglement spectrum of 
this system by cutting the cylinder at a; = 0. Thus the 

(n) 

projections ^ of the Landau level wavefunctions onto 
the right half-cylinder will also used to evaluate the mo¬ 
mentum polarization. We list their analytic forms here 
for the first three Landau levels: 


p^(o) _ 


ierfc 


kiB 


( 1 ) _ 1 -(kiE 


= 


^( 2 ) _ 


' fk£B^ 

\\lk£B 

\ a / 

) ^ 2 a 




^-{kiB/a) 


, ^ r f k(-B 

+ -erfc 


[a, a"*"] = 1 


a 



























6 


Now, the range of filled k states is determined by 
the length of the system; the last orbital is centered at 
X = ±KPg. We find that the viscosity derived from the 
sum over k converges to its expected continuum value 
when Lx > i.e. when each half of the cylinder 

is wider than a single wavefunction. We show the result 
of the viscosity calculation when successively filling up 
to the first three Landau levels in Fig. 1. We see that 
the Hall viscosity contribution from each Landau level 
converges to the established result® 


(n) 

Vh = 


n 1 

27r£| 2 



(16) 


which is the Hall viscosity contribution coming from the 
n-th Landau level. The convergence criterion is unsur¬ 
prising given our treatment of the cylinder’s edges. The 
edges are not sharp, rather the edge of the cylinder is 
some region defined by the width of the last occupied 
wavefunction. If the cylinder is narrower than the width 
of its edges then it is no surprise that the result does not 
converge properly. 

Using the correlation functions above, and Eq. (7), 
nearly the same results are obtained by the momentum 
polarization method with 


Ares = H 2 ^ 


n,fc 


Aoes = 1111' 


,ikAy 


n k>0 


The product over n spans the occupied Landau levels. 
The Hall viscosity can be calculated with Eq. (11) by 
computing Ares and Aqes at several values of Ly and ex¬ 
tracting the quadratic fit coefficient. The result is shown 
in the lower panel of Eig. 1, which shows that the calcu¬ 
lation converges when Lx > This is the same 

criterion as for the convergence of the momentum trans¬ 
port calculation: each half of the cylinder must be wider 
than a single wavefunction. 


IV. HOFSTADTER MODEL 

The utility of the two methods we have introduced is 
that they can be easily adapted to calculate the Hall vis¬ 
cosity in discrete, lattice systems in magnetic fields as 
well. The lattice systems have discrete translation and 
rotation symmetries, and have an additional length scale 
a, the lattice constant. Since there is not continuous rota¬ 
tion symmetry, then we can no longer appeal to the result 
that the viscosity is quantized in terms of the density. 
Furthermore, when considering momentum transport, we 
must consider the fact that continuous translation sym¬ 
metry is broken, and thus we are really considering the 
transport of quasi-momentum. Additionally, for the mo¬ 
mentum polarization technique, there is now a minimal 
Ay, i.e. the lattice constant in the y-direction. In lattice 



k / n 


FIG. 2. The spectrum of the Hofstadter Hamiltonian in 
Eq. (17) with q = 20. The first three Landau levels are high¬ 
lighted to illustrate the level filling scheme. 


systems we thus might expect that there is a maximum 
viscosity bound that is physically meaningful, i.e. when 
twisting the lattice by a single-lattice constant causes the 
transport of a full reciprocal lattice vector of momentum, 
then it is as if nothing has been transported. We will 
save a careful discussion of some of these issues to fu¬ 
ture work. For now we will compare the results of the 
two methods to see if they give matching results for the 
lattice viscosities, and moreover, if they both converge 
to the continuum limit when the magnetic length be¬ 
comes much longer than the lattice scale. As an aside, 
we note that Ref. 26 has also performed some viscos¬ 
ity calculations for the Hofstadter problem in a different 
context/methodology and recovers the continuum limit 
of the viscosity for small magnetic fields. 

We begin with the Hofstadter model,which is the 
tight-binding version of the integer quantum Hall prob¬ 
lem. The square lattice tight-binding model with rational 
flux (j) = p/q per plaquette has a Hamiltonian 

H = —tx(^n+l,ky^n,ky—'ty cos {ky — 2T:(j)n)(^^ j,^Cn,kyLLl.C. 

n, ky 

on a cylinder, where Cn^ky annihilates an electron in the 
y-momentum mode with wavenumber ky on the n-th site 
in the cc-direction. 

Although the Hofstadter model in the Landau gauge 
does not retain the fundamental translation symmetry of 
the lattice in the x direction, it is symmetric under trans¬ 
lation by a whole magnetic cell {q unit cells). To have pe¬ 
riodic boundary conditions in a torus geometry, we must 
respect this symmetry by having an integer number of 
magnetic cells, Nx = Iq for integer 1. Constructing the 
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FIG. 3. The Hall viscosity of the Hofstadter model calcu¬ 
lated by Eq. (3), with Ny = 51 and Nx = 2q — 1. The Hall 
viscosity of the continuum model is shown by a dashed line 
for comparison. 


cylinder from the properly periodic torus-by setting the 
wavefunction to zero on all the sites at one x coordinate- 
preserves the symmetry. Thus, the cylinder has the same 
number of unit cells as the torus, Iq, but it has one fewer 
site, Nx = Iq— 1."^® It is important to note that the con¬ 
sideration of the magnetic translation symmetry not only 
affects the construction of the lattice, but also the correct 
scaling of the viscosity (by the lowest Landau level den¬ 
sity). We also recall that with commensurate boundary 
conditions, the spectrum of the Hofstadter Hamiltonian 
on a cylinder has q nearly flat bands (Landau levels) con¬ 
sisting of Z — 1 states for each momentum mode ky. In 
the gaps between bands are edge states (one per mode 
ky) which connect the flat bands. We show an example 
of the energy spectrum with open boundary conditions 
and 9 = 20 in Fig. 2. 

Since we are only deforming the diagonal components 
of the metric, we can the standard nearest-neighbor tight- 
binding model given above. Hence, the lattice is de¬ 
formed through the hopping parameters and ty. Ab¬ 
sent the magnetic field, the tight-binding model bands 
are well-known: 

^(.kxi ky) - 2ta, cos (kx^x) cos )kyCly) . 

To see how the metric will enter, we can match parame¬ 
ters to the continuum Hamiltonian through a long wave¬ 
length expansion around the bottom of these bands: 

e(kx, ky) = ^ + ^(fcjOi)^^ . 


We can compare this to the Schrodinger equation: 

; ill 

‘ = + 

and then equate the coefficients of ki. Thus, we find the 
hopping amplitude in each direction is inversely propor¬ 
tional to the lattice constant: 

_ 

* 2maf 

Therefore, under the metric deformation in Eq. (2), we 
have tx cc c? and ty cx a~'^. Hence, we consider a de¬ 
formed Hamiltonian: 

^ ~ ^ cli+i^ky^n,ky 

n, ky 

- ^ cos {ky - 2TT(j)n)cl f^^Cn,ky + h.c.^ . 

(17) 

Now let us explicitly detail how the momentum- 
transport is calculated. The projected 

(quasi-)momentum is 

ul 

hky (m, ky\VR |m, ky) (18) 

ky m—1 

where Vr projects onto the right half of the cylinder: 

■Pi? = E I®) ■ 

31 — 0 

The integers m run over energy eigenstates at a given 
ky from 1 (lowest energy) to a value depending on the 
filling. As a reminder, we point out that if the v > q/2, 
the edge states associated with each Landau level above 
the middle of the spectrum are actually below the flat 
Landau level, rather than above, so one would need to 
be careful when choosing which states are filled if the 
viscosity of those Landau levels is of interest. The filling 
scheme for the first few Landau levels is illustrated in 
Fig. 2. We also note that, for our calculations, the site 
at which the cylinder is cut should fall on the boundary 
between magnetic cells, i.e. on a site n = rq for r S Z. 
We need to impose this condition so that the subsystems 
have commensurate boundary conditions, thus ensuring 
that the edge states of each half cylinder are the same as 
the physical edge states of the whole system.^® 

The viscosity can be computed directly from the eigen¬ 
states of the Hamiltonian (17) using Eqs. (3) and (18). 
At large q, the magnetic field is weak, and the magnetic 
length is much larger than the lattice spacing. In this 
regime we expect that the Hall viscosity should approach 
the continuum model. In fact, as in Fig. 3, we see that 
it does converge to the continuum result for the fillings 











we tested. As one increases the magnetic field, the ef¬ 
fects of the lattice will become more prominent. This 
figure also indicates that lattice effects more strongly af¬ 
fect higher Landau levels since the convergence to the 
continuum limit is slower. Eventually, as the magnetic 
field strengthens, i.e. as g —?> 0, the viscosity begins to 
depend on the lattice scale. From our results in the con¬ 
tinuum we expect that, when divided by the density, the 
viscosity should be a constant, independent of q. Instead 
we find that the viscosity has contributions that depend 
on q: 


2Tr£ 

h 


i-d' 


2^4 „( 2 ) 

h 


27r4 (3) 

h 


0.2499 

0.2499 

1.0042 

1.0042 

2.2289 

2.2289 


0.0017 0.3865 

0.0017 a 0.3865 

(-B 27r 

0.1513 4.3204 

0.1513 a 4.3204 

(-B 27r 

0.5938 2.2256 

^/q ^ q 

0.5938 a 2.2256 

(-B 27r 


where we have rewritten the q dependence in terms of 
the relevant length scales using the fact that qa^ = 2Tr£'^. 
We find that the viscosity is unchanged under B ^ — B, 
including the g-dependent contributions. 

The calculations in Fig. 3 are performed using deriva¬ 
tives with respect to the metric deformation a, but eval¬ 
uated at the isotropic point 0 = 1. For comparison, in 
Fig. 4 we fix g large enough (g = 120,180) so that the 
viscosities for the first three Landau levels have (nearly) 
saturated at the continuum limit, and then evaluate the 
momentum transport at different values of a. That is, we 
see how deforming around an initially anisotropic system 
affects the calculation. For values of a ^ 1 the sys¬ 
tem only has 180°-rotation symmetry, and is quasi-lD 
for large deviations. Fig. 4 shows that the viscosity of 
the Hofstadter model varies as a function of a itself; in 
comparison, the viscosity of the continuum Landau level 
is constant as a is varied. If the system is anisotropic 
we would expect the Hall viscosity to be controlled by 
more than one coefficient, for example ^ or 

^ 0- Helpfully, because of the a —>■ a~^ symmetry 
of the metric when we switch x and j/, we should be able 
to read off both viscosity coefficients from the same figure 
if we consider both a and a~^ simultaneously. However, 
we do not expect to enter the momentum trans¬ 

port calculation, and so we cannot extract that coeffi¬ 
cient from this figure. Finally, we note that if we directly 
compare the results at g = 120 and g = 180 as shown 
in Fig. 4, we find that dependence of the Hall viscosity 
on the anisotropy is weaker for larger g, which might be 
expected since lattice effects will naturally be less impor¬ 
tant when £b ^ a. In future work, it would be interesting 



FIG. 4. The Hall viscosity of the Hofstadter model at g = 120 
(red) and q = 180 (blue) as a function of a in comparison to 
the continuum quantum Hall model (dashed gray) for filling 
= 1, 2, 3. The system is isotropic when a = 1 and the unit 
cells are elongated in the y direction when a > 1. 


to see if any remnants of lattice anisotropy might survive 
to affect the thermodynamic limit. 

To close this section, let us compare these results with 
those obtained from the momentum polarization method. 
Note that in Eq. (18), the factor 

|w, ky) 

is just the aforementioned correlation function. From 
Eq. (7), we can use the correlation function to compute 
the momentum polarization phases: 

Ares= n \ [{l + e^^y^y) 

m,ky 

+ (1 _ (2 (m, kylVB. \m, ky) - 1)] 

Aoes= n \ + 

m,ky 

-h (l - (20 ((to, ky\X \m, ky)) - 1)] , 

(19) 

where 0 is the Heaviside step function, and x is the x- 
coordinate operator. On a lattice. Ay must be an integer 
in units of the lattice constant. The resulting viscosity 
calculation is shown in Fig. 5. The Hall viscosity ob¬ 
tained for the first three Landau levels agrees with the 
continuum value in the weak field limit. Where the calcu¬ 
lation converges, i.e. g > 20, it agrees qualitatively with 
the momentum transport method, although the momen¬ 
tum polarization calculation appears to deviate less from 
the continuum Hall viscosity at small g. 
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FIG. 5. The Hall viscosity of the Hofstadter model calculated 
by Eqs. (11), (19), with = 2q—l. The Hall viscosity of the 
continuum model is shown by a dashed line for comparison. 


V. CONTINUUM DIRAC LANDAU LEVELS 


For our second set of examples we focus on 2D Dirac 
fermions in a magnetic field. The quantum Hall effect 
in this type of system became fundamentally important 
with the rise of graphene,and more recently has be¬ 
come relevant in the study of 3D topological insulators 
with low-energy surface fermions of Dirac nature.^® We 
can describe the Landau level problem in the (massive) 
Dirac Hamiltonian under shear strain with the Hamilto¬ 
nian 


H = a {px + eA^) + a ^ {py + eAy) + ma^ ( 20 ) 

where we have again chosen the a-dependent met¬ 
ric/frame in Eq. (2), and cr“ are the usual Pauli matrices. 
As above for the Schrodinger equation, let us consider a 
cylinder which is periodic in the y-direction with circum¬ 
ference Ly, and in the Landau gauge where A^ = 0 and 
Ay = Bx for a uniform magnetic field normal to the sur¬ 
face. Again the Hamiltonian can be written in terms of 
the raising and lowering operators ( 12 ): 


H = 


m '^2heBa \ 
y/2heBa^ —m ) 


The Landau level wavefunctions are, up to normalization. 






v^AT‘ J 


Hh/^Po 
1.2 - 


1.0 


0.8 - 


0.2 


n=±2 


: • 


.it* 


4 


Lx 


Ib 


FIG. 6. The Hall viscosity (ph) of the specified integer 
quantum Hall states of the Dirac model (20) calculated by the 
momentum transport method. As in Fig. 1, the calculation 
converges when > 2y/^£B- The derivatives in Eq. (3) 
were taken numerically with (PyVR) given by Eq. (21). The 
Hall viscosity is given in units of hpo where po = ll2n£% is 
the electron density of the lowest Landau level. 


where 7 = m/V^heB is the ratio of the two energy scales 
in the problem, ^ are the Schrodinger Landau level 
wavefunctions wavefunctions (13), Py = hk, and we have 
denoted for convenience p±ni'l) = 7 ± 1 / 7 ^ -I- n. The 
energies of each Landau level are 

Eq = —jV2HeB = —m 
E±n = iV 2heB \/ 7 ^ + n. 


Now that we have the Landau-level wavefunctions, we 
can calculate the viscosity using Eq. (3), i.e. by differen¬ 
tiating 


K 




(0) 


{Py^R. 


(n^O) 


k=-K 
K 

= E 

k^-K 


+P±n(7)^C'E 

n + p±n{lf 


( 21 ) 


with respect to a^. Recall that defined in Eq. (15), 

and matches our earlier results since the Dirac Landau- 
levels are constructed from the Schrodinger Landau- 
levels. 

Because of the connection between the Dirac and 
Schrodinger Landau-levels we conclude that the Hall vis¬ 
cosity of each Landau level in the continuum Dirac sys¬ 
tem is given by 


„( 0 ) 

^H,S 


n = 0 



n 0 

n+p±. 



( 22 ) 


(n) 

where g is the Hall viscosity of the nth Landau level of 
the continuum Schrodinger equation as given in Eq. (16). 
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In the massless limit, when m = 7 = 0 so that p±n = 
±\/n, we find 


(n) 

'^H,D 


h/ n = 0 

h \n\ / ( 47 r^^) n ^ 0 . 


(23) 


This result is in agreement with previous work by 
Kimura^® based on an adiabatic curvature calculation, 
except for the n = 0 level for which we have found a 
value twice as large. We attribute the difference to a 
probable error in the normalization of the zeroth Lan¬ 
dau level in Ref. 15. We confirm the results numerically 
in Fig. 6 using the momentum transport method of cal¬ 
culating the Hall viscosity. Because the Hall viscosity 
of the Dirac Landau levels is expressed in terms of the 
Hall viscosity of the Schrodinger Landau levels, the con¬ 
vergence criterion is expected to be the same. Indeed, 
we find the result converges to the expected value when 
Lx > The same result is obtained by the mo¬ 

mentum polarization method, with similar convergence 
criteria, though we do not show the figure here. 

Let us now test if the Dirac calculation reproduces the 
Schrodinger result in the large mass limit. Thus, we will 
consider the 7 —>■ ±00 limit. In either limit, 

, n 

P+„(7)«27+ — 

and the resulting wavefunctions are 


^-(n) ^ ^ , 


The limiting values of the wavefunctions can easily be 
determined by considering the order, with respect to 
7 , of each component of the spinors. Additionally, the 
wavefunction is completely unmodified in this limit. 
From this result we can conclude immediately that, in 
the infinite mass limit, the Dirac Landau levels carry the 
same set of values of the viscosity as Schrodinger Landau 
levels, though we still need to see how they are organized. 
Additionally, in this limit the energy eigenvalues are 

Eq = —m 

E±n « EV^heB ^7 -I- —^ = ilAol ± huin 

with oj = \eB/m\ the usual cyclotron frequency. The 
spectrum has a gap of width 2|Ao| with Landau lev¬ 
els above and below separated from neighboring Lan¬ 
dau levels by a gaps of uniform width hui, much like the 
Schrodinger spectrum. 

The conclusions so far hold generically in the 7 —> ±00 
limits. Let us now consider each limit independently, and 


m = 


0 


m = 0 






FIG. 7. The spectrum of the lattice Dirac Hamiltonian 
in Eq. (24) at m = 0 (top) and m = 4 (bottom) The n = 
0,±1,±2 Landau levels are indicated. A zoomed-in view of 
these Landau levels is presented in the right panels of each 
row. 


furthermore, let us consider taking each limit by fixing 
B and sending m —t ±00, respectively. In either case, 
the wavefunction of the n = 0 Landau level is essentially 
unchanged from the Schrodinger system. When m —> 00, 
the n = 0 Landau level sits at the top of the valence 
{E < 0) band, separated from the n > 0 Landau levels 
by the (large) mass gap. On the other hand, when m —>■ 
—00, the n = 0 band sits at the bottom of the conduction 
{E > 0) band with only the cyclotron gap separating it 
from the n > 0 states. It is this configuration, when 
m —t —00, and with the E < 0 states filled, which more 
precisely matches the Schrodinger case. This should not 
be surprising; the ma^ term of the Hamiltonian attaches 
a positive mass to the n = 0 Landau level when m < 0. 
Thus, we see that the massive Dirac case matches the 
Schrodinger case if one focuses on the positive energy 
levels when to —> — 00. 

Now that we have discussed some properties of the 
continuum Dirac model, let us consider a lattice version. 


VI. LATTICE DIRAC MODEL 


Let us consider a lattice regularization of the con¬ 
tinuum Dirac model. Despite the fact that the Dirac 
Landau-level spectrum is celebrated because of its appli¬ 
cation in graphene, we will not consider such a honey¬ 
comb lattice model. The reason is that they honeycomb 
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model presents extra difRculties. For example, there are 
not only multiple Dirac cones, but the cones are located 
away from the F-point in the Brillouin zone. The lat¬ 
ter issue leads to results which are not easily compara¬ 
ble with the Dirac viscosity calculation in the continuum 
limit. We have performed cursory calculations on such a 
system, but we will leave the discussion of lattice viscosity 
calculations when the low-energy states are near generic 
points in the Brillouin zone to future work. Instead we 
will consider a simpler lattice model for a Dirac fermion 
on a square lattice. When the metric deformation of 
Eq. (2) is included, the lattice Dirac model Hamiltonian 
(on a square lattice in a cylinder geometry periodic in 
the y direction, with rational flux (j) = p/q per plaquette) 
is 

-bh.c.) 

n,ky 

+ 'ii,ky sin {ky - 2TT(l)n)a~^a^Cn,ky 

+ 4ik [2 - m - cos {ky - 27r(/)n)] , 

(24) 

where Cn,ky is a two-component annihilation operator. 
This model has a single gapless Dirac cone when m = 0 
or TO = 4. For to = 0 (to = 4) the Dirac cone is located 
near k = (0,0) (k = (7r,7r)) in the Brillouin zone. Like 
the Hofstadter model, since we have included a magnetic 
field, the system geometry is chosen to preserve the mag¬ 
netic translation symmetry with an integer number of 
magnetic cells that are q sites wide in the x direction. 
Since we are using a cylinder geometry the lattice should 
have Nx = Iq — 1 so that the boundaries are commensu¬ 
rate.^® The Landau levels at to = 0 and to = 4 are shown 
in Fig 7. We immediately recognize the similarities in the 
two cases, but should point out a major difference, i.e. 
that the edge states for to = 0 (to = 4) are located near 
ky — 0 {ky = tt). Below we will discuss how this differ¬ 
ence affects the results for these two cases. 

The viscosity can again be calculated by projecting the 
total momentum onto the right half of the cylinder, as in 
Eq. (3), and then differentiating with respect to c?. The 
momentum projection is 

{PyVi:^ = ^^kky (j, ky\VB. li, ky) 

ky jocc. 

where Vr projects onto the right half of the cylinder: 

Nj: 

2 

'Pr = '^ |x,cr) (x,cr|. 

x=0 cr=±i 

The integers j run over the occupied energy eigenstates 
at a given ky. For most of the cases we consider we only 
fill the Landau levels near half-filling. We note that near 
half-filling the n-th Landau level consists of the states 



FIG. 8. The Hall viscosity of the lattice Dirac Hamilto¬ 
nian (24), calculated by Eq. (3) with Ny = 51 and Ab, = 2q—l. 
The indicated individual Landau levels are filled. The dotted 
grey line indicates the Hall viscosity of the continuum Dirac 
Landau level given by Eq. (22). 



FIG. 9. The Hall viscosity of the lattice Dirac Hamilto¬ 
nian (24), calculated by Eq. (3) with Ny = 51 and AL, = 2q—1. 
The Landau levels are filled from the bottom of the spectrum 
through the indicated level. Because filling this way causes 
the number of filled Landau levels to vary with q, a linear 
term O.Ollg has been subtracted from each series. 
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FIG. 10. The Hall viscosity of the lattice Dirac model Hamil¬ 
tonian (24), calculated by Equation (3) with Ny — 51 and 
Nx = 2q—l. The n = 0 Landau level Hall viscosity is plotted 
in each case. The m = 4 Hall viscosity is plotted twice. The 
unshifted plot shows the viscosity when the Brillouin zone is 
unshifted so that the Dirac point is at fc = ±7r. In the shifted 
plot, the Brillouin zone has been shifted so that the Dirac 
point is once again at fc = 0, in which case the Hall viscosity 
agrees exactly with m = 0. 


I I I i I I 

* -W -M- *- *- *- *-*-*-) 


_ X 


X m = 0 • 

X m = 4 (unshifted) x 
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X 


X 
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(see Fig. 7) 


j G 


{^x + nl: Nx + (tJ + 1)Z] 
{Nx + (n - 1)1, Nx + nl] 


m = 0 
TO = 4. 


Notice that the 0-th Landau level moves from the bottom 
of the conduction band at to = 0 to the top of the valence 
band at to = 4. This is clearly shown in Fig. 7. 

The Hall viscosity of the lattice Dirac model was cal¬ 
culated at TO = 0 by the momentum transport method 
in Eq. (3) to obtain the results in Fig. 8. The values 
here represent the viscosity calculations from individu¬ 
ally filling (not successively filling) the n = 0,±1, and 
±2 Landau levels (where n = 0 is referenced to the ze¬ 
roth Landau level of the Dirac point, not the bottom of 
the entire bandwidth). To help illustrate, we have shown 
which Landau levels were filled in Fig. 7. We note that 
the lattice calculation converges to the continuum value 
in the large-g (weak magnetic field) limit, i.e. the Hall 
viscosity of the lattice system approaches the continuum 
value in the limit where the magnetic length £b is much 
larger than the spacing between unit cells. As the mag¬ 
netic field strength increases, so does the effect of the 


lattice, with the viscosity taking on ^-dependent terms: 


27r4 (_ 2 ) 
h 

- 0.9868 -b 

27^4 „(-i) ^ 

- 0.5018-b 

(0) _ 
h 

- 0.2498-b 

27r4 (1) 
h 

- 0.5025-b 

2’^^! ( 2 ) 
h 

- 0.9802 -b 


0.4276 a 12.0267 

Ib 2tt 

0.0576 a 1.7067 

1 /^ £b 27r £g 

0.0045 a 0.8290 

A/27r £b 27r 

0.0857 a 1.3075 

£b 27r 
0.6299 a 14.1092 

£b 27r £‘^g’ 


where we have used the relation = 27r^^. As pre¬ 
dicted by the continuum calculation, the Hall viscosity 
converges to approximately the same value for positive 
and negative Landau levels. 

Now let us consider what happens if we fill the Landau 
levels from the absolute bottom of the spectrum, instead 
of filling individual Landau levels near the Dirac point 
as was done in Fig. 8. We show the results in Fig. 9. 
Interestingly, we found that the momentum transported 
scaled linearly in q, and we subtracted off this contri¬ 
bution to make Fig. 9. This linear scaling should be 
expected because the number of filled Landau levels is 
proportional to q when filling from the bottom of the 
spectrum. After subtracting off the linear term O.Ollg, 
as shown in Fig. 9, we find that the momentum trans¬ 
ported indeed saturates to a fixed value in the low-field 
limit. However, the saturation values do not match the 
continuum results. Instead, it is only the differences in 
the viscosities between filling the n-th and (n -I- l)-th 
Landau levels that exactly matches the continuum value 
for the viscosity of the added Landau level. This is a 
surprising result, as it indicates that for lattice systems 
in magnetic fields the magnitude of the viscosity may be 
somewhat regularization dependent, but the difference in 
viscosities seems to retain a more universal character. It 
would be interesting to see if this is a generic feature, an 
artifact of this model, or can be attributed to finite size 
effects. 

Let us now consider the other massless limit of this 
model when to = 4. We show the viscosity, calculated 
via momentum transport, in Fig. 10. The bare result for 
the viscosity shows a monotonically decreasing function 
that does not converge for large q. However, to properly 
interpret this result, care must be taken to recenter the 
Brillouin zone. If we keep to fixed, but send fcy —)■ fcy — tt 
in the Hamiltonian, then the momentum transport cal¬ 
culation exactly recovers the result at to = 0. If the Bril¬ 
louin zone is not shifted, extra momentum is transported 
since the edge states are located near ky = tt, and this 
leads to a different result for the viscosity. We show the 
TO = 0 case, as well as the shifted and unshifted results 
for TO = 4, in Fig. 10. Generically, when the low-energy 
Dirac point (s) are away from k = 0 in the Brillouin zone 
there is extra momentum transport due to the overall 
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FIG. 11. The Hall viscosity of the lattice Dirac Hamilto¬ 
nian (24), calculated by the momentum polarization method 
with Nx = 2q—l. The indicated individual Landau levels are 
filled. Where points are missing it indicates a failure of the 
fitting required to extract the viscosity. 


momentum shift of the cone. While we have been able 
to adjust the calculation for this simple case (and conse¬ 
quently any case with a single Dirac cone), the question 
of how to compare the viscosity of multiple Dirac points 
at generic momenta to the continuum limit remains a 
topic for future work. 

The Hall viscosity of the lattice Dirac model was also 
calculated by the momentum polarization method as in 
Equation (8). As in the Hofstadter model, the correlation 
function is 


= ky\T^R\j, ky). 

Equation (7) allows us to compute Ares and Aqes, 

Ares = H ^ 

+ (2(j, kylVnlj, ky) - 1)] 

Ages = H ^ 

J,ky 

+ (20((j, ky\x\j, ky))-l)]. 

Again, 0 is the Heaviside step function, and x is the 

x-coordinate operator. As with any lattice model. Ay 
must be an integer in units of the lattice constant. The 
Hall viscosity obtained this way agrees with the momen¬ 
tum transport calculation, showing the same convergence 
to the continuum value of rj^ at large q (Figure 11). 
Note that although both methods show a deviation from 


the continuum Hall viscosity at small q, the momen¬ 
tum polarization method shows a smaller deviation and 
with opposite sign. Points are missing from these figures 
where the fitting required for the momentum polarization 
method has failed. 


VII. DISCUSSION AND CONCLUSION 

We have applied two techniques for calculating the Hall 
viscosity in integer quantum Hall systems. Both meth¬ 
ods seem to capture similar results for continuum models, 
and, more interestingly, both were successfully applied 
to lattice models. Our original momentum transport 
method gives results in agreement with the momentum 
polarization method previously described, especially in 
the weak magnetic field limit. While we have seen that 
there are lattice-scale dependent contributions to both 
the momentum transport and the momentum polariza¬ 
tion at strong magnetic field, these corrections seem to 
be method dependent, at least for the system sizes and 
parameters we have chosen. 

We have demonstrated that either method can deter¬ 
mine the Hall viscosity of an isotropic system. However, 
further work will be required to fully characterize the Hall 
viscosity of anisotropic systems, though we demonstrated 
that two of the three non-vanishing Hall viscosity coeffi¬ 
cients can already be computed using our methods and 
geometry. We found that the Hall viscosity coefficients 
were dependent on the amount of anisotropy, though it is 
unknown if this variation survives in the thermodynamic 
limit. 

Finally, we have calculated the viscosity for the Dirac 
Landau-level system, and shown that there is a relation¬ 
ship between the Hall viscosity of continuum Schrodinger 
and Dirac Landau levels. Namely, that the latter recovers 
the viscosity of the former in the infinite mass (m —> oo) 
limit. Our methods reproduce the known results for these 
models. Furthermore, we show that both the Hofstadter 
and lattice Dirac models approach the appropriate con¬ 
tinuum Hall viscosity as the magnetic field H —>■ 0, with 
deviations at stronger fields that depend on the lattice 
scale. Our results for the lattice Dirac model suggest 
that it may be the viscosity difference between Landau 
level fillings which is actually quantized (in units of den¬ 
sity) in lattice regularized models. We uncovered some 
difficulties in treating lattice Dirac systems with multiple 
Dirac points, when they are located at generic points in 
the Brillouin zone; further work will be required to treat 
some systems of interest, such as graphene. It is also of 
interest to understand the competition between time re¬ 
versal breaking arising from the applied magnetic field, 
and an intrinsic time reversal breaking coming from a 
massive Dirac model, i.e. a Chern insulator. The latter is 
also expected to have a non-vanishing, field-independent, 
contribution to the viscositythough it has yet to be 
calculated in a lattice regularization. The article Ref. 49 
addresses some aspects of this last topic. 
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